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\q We discuss models that can account for today's dark energy. The underlying 

CN cosmological constant may be Planck scale but starts as a redundant cou- 

, , pling which can be eliminated by a field redefinition. The observed vacuum 

energy arises when the redundancy is explicitly broken, say by a non-minimal 
coupling to curvature. We give a recipe for constructing models, including 
R + 1/R type models, that realize this mechanism and satisfy all solar sys- 
tem constraints on gravity. A similar model, based on Gauss-Bonnet gravity, 
provides a technically natural explanation for dark energy and exhibits an 
interesting see-saw behavior: a large underlying cosmological constant gives 
\q rise to both low and high curvature solutions. Such models could be statis- 

ts tically favored in the string landscape. 
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1 Introduction 



The present-day universe is undergoing exponential expansion on a timescale 
set by today's Hubble parameter H ~ lCT 33 eV. One of the biggest chal- 
lenges in physics is understanding why Ho is so small relative to the energy 
scales of particle physics. This is the puzzle of the cosmological constant, 
reviewed in [H [21 El 13 El- 
Many solutions have been proposed [Bj. An approach which is currently 
popular is to invoke landscape / anthropic ideas: in a potential with enough 
minima, at least one should have the observed vacuum energy [7j. Although 
this may ultimately prove to be the right explanation, it's still important to 
explore alternatives. In fact, the existence of a landscape may require us to 
explore alternatives. To explain this point, let us note that in a simple model 
of flux compactification one has the potential [8] 



Here — Vq is a large fixed negative energy density, N is the number of non- 
trivial cycles on some compactification manifold, the Cj are order one con- 
stants, and the integers measure the fluxes through the various cycles. 
The number of flux vacua then grows with the vacuum energy roughly as 



For large N the number of vacua grows extremely rapidly with energy. Now 
suppose there is a mechanism for (nearly) cancelling a large positive vacuum 
energy, leaving a small effective cosmological constant. We do not expect 
to find such a mechanism without introducing small parameters [I]. But in 
view of the growth (IT]), such a mechanism could be statistically favored in 
the landscape, even though it requires carefully adjusted couplings. 

This leads us to consider theories which nearly cancel a large underlying 
vacuum energy. The basic mechanism we shall explore is the following. Sup- 
pose there is a potential V for a collection of fields ip. We assume that V(ip) 
vanishes somewhere in field space, but (as should be generic) we assume the 
derivatives of V are non-zero at that point: 



N 



i=l 




(1) 
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dV 



7^ for some i. 
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Then in a neighborhood of ipo we can redefine fields, introducing a new field 
4> such that 

V = Pvac -<j). (2) 

Here p vac is a constant representing a large underlying vacuum energy. Clearly 
at the level of the potential (J2J), p vac has no meaning: it is a redundant cou- 
pling which can be eliminated by a field redefinition 

-> + a p vac -> p vac + a . (3) 

However generically this pseudo-symmetry will be explicitly broken by other 
terms in the Lagrangian. By choosing appropriate breaking terms, we will 
be able to write down models that have de Sitter solutions of arbitrarily 
small curvature, even in the presence of an underlying Planck-scale vacuum 
energy. The necessary breaking terms will involve extremely small couplings. 
In practice this allows us to shuffle fine tunings away from the cosmological 
constant and into other sectors of the theory. Such tunings may be expected 
to arise in the landscape, as we emphasized above. Also by shifting small 
parameters around, we will be able to construct models in which a small cos- 
mological constant is technically natural (stable under radiative corrections). 
Interestingly, these models also exhibit a kind of see-saw mechanism, in that 
a large underlying cosmological constant gives rise to multiple solutions, some 
with high curvature, some with low curvature. 

An outline of this paper is as follows. In section [2] we develop the idea of 
making vacuum energy nearly redundant. In section [3] we present a model- 
building recipe based on fields with non-minimal couplings to the scalar cur- 
vature. In section [4] we discuss a particular example in more detail: a model 
with an Rcf 2 coupling that can be thought of as a generalization of R + 1/R 
gravity j9]. In section [5] we discuss another example, based on Gauss-Bonnet 
gravity, which provides a technically natural explanation for dark energy. In 
the appendices we show that all models we consider can be made compatible 
with cosmological and solar system tests of gravity, and give further details 
on issues of classical and radiative stability. Readers who are interested in 
the concrete examples could skip directly to sections [4] or [5j which can be 
read independently of the rest of the paper. 

Conventions: our metric signature is ( — h++). Quantities which appear 
in the underlying Lagrangian include the bare gravitational coupling k = 
1/mp and the bare vacuum energy p vac . Quantities measured at low energies 
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include the effective gravitational coupling n eS = l/rrfp = (SttGn) 1 ^ 2 ~ 
10 18 GeV and the effective vacuum energy p^ c ~ (lCT 3 eV) 4 . 



2 Redundant vacuum energy 

Conventionally one imagines that vacuum energy arises from an (effective) 
potential V(ijj) for a collection of matter fields ipi. To find static solutions 
one goes to a minimum of the potential V'(ip) = and finds the vacuum 
energy 

Pvac = V(l/>)\y l=0 . 

Note that p vac is invariant under field redefinitions of the ipi. The effective 
action is 




and one has the usual tuning problem for vacuum energy: p va _ c / m p "C 1- 

As an alternative to the conventional scenario imagine a potential which 
crosses through zero somewhere in field space. That is, suppose for some ipo 

V(iJo)=0 but V'(^ )^0. 

We can redefine fields in a neighborhood of ipo to set V((f>) = p vac — J (p. (The 
potential may have critical points elsewhere in field space, in which case this 
redefinition can only be made locally.) For simplicity we assume a canonical 
kinetic term for - an assumption we will relax in the next section - and 
consider the action 

S = j d A x^g- (J^R - X -d^4> + J<P- p vac ^) . (4) 

The interesting thing about this action is that it's invariant under 

-» + a p vac -> p vac + Ja. 

This means the bare vacuum energy p vac has no observable consequences: 
it's a redundant coupling which can be eliminated by a field redefinition [TP] . 
This is unlike the source J which could be absorbed by rescaling the field 
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but would then reappear in the kinetic term. Of course Q has no static 
classical solutions, rather there are only cosmological solutions in which the 
scalar field rolls down the hill. p vac can be absorbed into initial conditions 
for the field - say the value of the field at which 8q4> vanishes. 

To obtain de Sitter solutions one could imagine deforming Q by turning 
on other couplings. For example one could consider the theory 

S = Jd 4 x^9 (f(<l>)R - Id^dy - V{<j>)\ . (5) 

Now we can find de Sitter solutions. Indeed for constant fields the scalar 
equation of motion requires 

dj R = ^ he - R = w (6) 

and the Friedmann equation is 

R = TV If . (7) 

Provided dV/ dip and df / d<p are non-zero these equations have novel de Sitter 
solutions: novel in the sense that (unlike the conventional scenario) the effec- 
tive vacuum energy is not determined by critical points of the potential, and 
(unlike quintessence) our models do not invoke a slowly-rolling scalar field. 
We will refer to vacuum energy in our models as being "pseudo-redundant." 

What does making vacuum energy pseudo-redundant achieve? We should 
state at the outset that it does not solve the fine-tuning problems associated 
with the cosmological constant. As we will discuss in more detail in the next 
section, the model ([5J will have an effective Planck mass m e p and an effective 
vacuum energy p°f c , and we require 



(mf) 2 (mfy 

for the model to be compatible with observation. This can be achieved, but 
only by tuning parameters, in accord with Weinberg's no-go theorem pQ. 
However, as is manifest from the construction, by making vacuum energy 
pseudo-redundant we have shifted the necessary tunings away from the un- 
derlying cosmological constant and into other sectors of the theory. This 
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freedom may be important in the landscape. Also, as we will see in section 
[5J it will allow us to build technically natural models for vacuum energy]]] 



3 An F(R) model-building recipe 



In this section we consider models which generalize ^ by including a non- 
minimal kinetic term for the scalar, 

. 1 



s = J d A xV=g U{(f>)R - - v{(j>)j + s mattcr (8) 

where matter is minimally-coupled to g^ v . Models of this type, with appli- 
cations to dark energy, have been reviewed in [Til 112] . We will refer to them 
as F(R) models^ 

As in the previous section, de Sitter solutions obey 

dV 

R = — and R = 2V/f . (9) 
ctj 

^From the action we can read off the effective Planck mass 

{mff = 2</(0)> (10) 
and the effective vacuum energy 

p c l = (vm. (11) 

We want to choose the functions / and V so that these quantities take on 
their observed values. However we also want to construct models in which the 
de Sitter solutions are (meta-) stable and compatible with solar system tests 



^ull disclosure: models such as (J5J) can be put in Einstein frame via a conformal 
transformation, see appendix |A.1| and in Einstein frame the necessary tuning has the 
conventional cosmological-constant form. This is further motivation for considering models 
which do not have an Einstein frame description, such as the Gauss-Bonnet model of 
section [5| 

2 In the literature F(R) gravity is usually taken to mean no scalar kinetic term, i.e. 
h = in As we discuss below, our models agree with F(R) gravity on cosmological 
solutions but differ inside the solar system. 
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of gravity. One way to achieve this is to make h large. This does not affect 
the de Sitter solutions, but it suppresses any spatial or temporal variation of 
4>, so for sufficiently large h the model passes all solar and stability tests and 
does not significantly alter the cosmological history during matter/radiation 
domination. Stability is studied in more detail in appendix |A.1[ solar tests 



are analyzed in appendix A. 2, and cosmology during matter/radiation dom- 



ination is studied in appendix A. 3 



^From now on we will assume h is large enough to pass solar system, 
cosmology, and stability tests, and we turn our attention to constructing 
models that yield the desired low curvature de Sitter expansion. We could 
analyze this using the scalar-tensor form ([8]). But it is convenient to assume 
constant <fi and integrate out the scalar field, to obtain the effective action 
for gravity 

S = Jd 4 x^C(R) C(R) = fR-V. (12) 

Here / and V are implicitly functions of the curvature scalar, determined by 
solving R = dV/df to find / = f(R)- This is simply a Legendre transform^] 
The transformation can be inverted with [131 E] 

f = dC/dR V = fR-C (13) 

to go back to the scalar-tensor description. As usual V here can be expressed 
purely in terms of / using the relation / = dC/dR. Strictly speaking this 
gives us the scalar-tensor action in terms of R and /; we can then imagine 
an arbitrary functional relationship between / and (ft. 

It's convenient to write the effective gravity action in terms of 

y = R/4(mff (14) 

and setJ3 

C(R) = {mffF{y) . (15) 

What conditions do we need to impose on F(y) to obtain a low-curvature 
de Sitter solution? The Friedmann equation ^ fixes R = 2V/f, so we must 
have 

( m fYF(y) = fR-V = V. 

3 Intcgrating out the scalar is analogous to starting from a first-order Lagrangian pq — H 
and integrating out the momentum using the equation of motion q = j^- . 



Here mfF is the fixed constant defined in (10 1 
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This equation should hold at y = e, where e = R /4:(m e p) 2 , R being the 
curvature today. Identifying (V) = p^ c = e(m e p) A we have the condition 
F(e) = e. This will give the right radius of curvature in Planck units. But to 
get the right effective Planck mass we also need (/) = \{rrfp) 2 which from 



the first equation in (13) means that F'(e) = 2. 

To summarize, any function F(y) with 

F(e) = e F'(e) = 2 (16) 

will produce a low-curvature de Sitter solution with the right effective Planck 



mass. Clearly there are an infinite number of functions which satisfy ( 16 ) . We 
present a number of examples in appendix [Bj and study a particular example 
- the R(f) 2 model - in more detail in the next section. Some of these examples 
might come as a surprise. For instance, models such as C ~ (1/-R) + (V-R) 2 , 
where an Einstein Hilbert term is completely absent, are actually viable 
both in the cosmological and the solar system contexts. Interested readers 
are urged to consult the appendix. 

Here we make some general comments on tuning issues. Any function 



satisfying (16) is fine-tuned, in that its value near zero is so small. This crit- 



icism applies, for example, to ordinary gravity with a cosmological constant, 



which corresponds to F(y) — 2y — e. This particular solution to (16) does 
not interest us so much, because we are more interested in models in which 
F(y) has an 0(1) constant piece representing a Planck-scale bare vacuum 
energy. 

F(y)~-1 + -.- (17) 
Satisfying both ( [16]) a nd ( [T7| ) requires an even greater degree of fine-tuning 



than usual, since (16) by itself suggests that F(0) « — e, while (17) by itself 
suggests that F(0) ~ —In That is, in F(R) gravity it is unnatural to have a 
small effective vacuum energy. But it is doubly unnatural to obtain a small 
effective vacuum energy by nearly cancelling off a large underlying vacuum 
energy: to achieve this requires an even greater degree of fine-tuning than is 
usually associated with the cosmological constant. Nonetheless we present 
some explicit examples of F(R) gravity with Planck-scale bare vacuum energy 
in appendix |Bj and we discuss one model in more detail in the next section. 
Readers who are concerned with tuning issues are advised to skip to the 
Gauss-Bonnet model of section |5l 

5 As in appendix [b] one could consider models with inverse powers of y. However this 
does not help with tuning as in such models F(0) diverges. 
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It is also worth commenting on the difference between our F(R) models 
and the conventional F(R) models often found in the literature (e.g. (9j [15] 
and references therein). The latter models set the bare vacuum energy to zero 
by hand (for instance p vac would be set to zero in the Rcf 2 example we consider 
next). Moreover these models have no kinetic term for the scalar field in 
Jordan frame, i.e. h = in ([8]). This implies that even for the same F(R), 
our predictions for subhorizon fluctuations would be quite different from the 
conventional model, even though we agree on the cosmological background. 
In general, our predictions are much closer to those of general relativity, 
although there could be detectable deviations in solar system tests. 

Finally, let us point out that the model-building recipe we have presented 
can be easily extended. For example, one natural extension is to demand that 
the solution be stable. In terms of the Einstein-frame potential V = ( 2 ( K g)2/)2 
defined in (145 ), stability can be achieved by having V" > 0. Evaluated on 
the solution (9j, V" > is equivalent to V" > 8e. Noting that V" = 
1/ (d 2 C/dR 2 ), the condition for stability is simply < F"(e) < 2/e. More 
generally, a program for reconstructing F(R) gravity from a given expansion 
history has been developed [121 HB] • 



4 Rcj) 2 model 

In this section we consider a particular example of F(R) theory, the R<p 2 
model with action 

S = j d A X^ (~R - h^4> - ^R4> 2 + J<f>~ Pvac) • (18) 

We derive this action from our general recipe in appendix |Bj but with the 
couplings k, £, p vac already adjusted so as to give the correct de Sitter radius. 
Here, rather than use our previous results, we will analyze the model from 
scratch. Note that we have set h, the kinetic coefficient in ([8]), to unity. A 
large h is equivalent to a small £ through a field redefinition. 

The action has a term linear in which involves a source J. We will 
not specify the dynamics that give rise to J, aside from noting that it could 
arise from a Yukawa coupling ipifxp with a condensate (tpip) 7^ 0. We have 



8 



included a non-minimal coupling to curvature |£i?0 2 . Also we have tuned 
the underlying <p mass to zero; as can be seen below, we require m| < £if 2 for 
the model to work. More generally any additional potential for <f), if present, 
would have to be extremely flat. Finally, particle physics could be taken 
into account by minimally-coupling the standard model to g^; any particle 
physics contributions to the vacuum energy can be lumped into p vac . 

To find cosmological solutions it's convenient to solve for <fi in terms of 
the source, <ft = J/£,R where we have assumed that <fi is constant. Plugging 
this back into the action we get 

Aside from the vacuum energy term, this is the generalized gravity model of 
Carroll et. al. [9]. The Friedmann equation for this model is 

1 k 2 T 2 

H 2 — - 
11 — „ rf Pvac 



where we have assumed a constant scalar curvature R = 12H 2 , appropriate 
for de Sitter space in a flat FRW slicing. This quadratic equation for H 2 has 
solutions 

„2 



H 2 = - 
6 



, 2 3J2 
Pvac±<l Pt 



vac 



(20) 



A£k 2 

Suppose J 2 ^ 2 Pvac- Then there is the conventional or high-curvature 
solution in which the Einstein-Hilbert term dominates and H 2 w |^ 2 p V ac- 
But there is also a low-curvature solution in which the 1/R term dominates 
and H 2 f» J 2 /16^p vac . For ^p va c > this provides a mechanism by which a 
large underlying vacuum energy can drive a slow Hubble expansion. We will 
refer to such pairs of high- and low-curvature solutions as see-saw solutions 
for dark energy]^] 

It is important to note that to obtain the effective R+l/R description of 
the R<p 2 model we had to integrate out a light scalar with a mass set by R. 



So the action (19) provides an effective description of cosmology but cannot 



6 Similar low-curvature solutions, driven by a large matter energy density, were found 
in |15] . In our general model-building recipe of section [3] there is no real need to demand 
that solutions come in see-saw pairs, although such pairs could arise. 
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be used on sub-Hubble distances, where one must revert to the underlying 
scalar-tensor theory (18). Indeed by plugging a constant vev for <p into the 
action ( 18 ) we can read off the effective Planck mass 

K ff ) 2 = ^-£<0> 2 (21) 

and the effective vacuum energy 

P 6 ! = Pvac - J{il>) ■ (22) 

Assuming constant <p, the Friedmann equation is 

^-£0 2 W = ^(Pvac- Jcfy (23) 



and the cf) equation of motion fixes (f) = J/12C,H 2 . Note that (23) can be 
rewritten as 

[mp ) ti - -p vac 

which is the usual Friedmann equation in the presence of vacuum energy. 
So in terms of the underlying R<p 2 theory, given m e p ~ 10 18 GeV and H ~ 
10~ 33 eV, the scalar field will adjust its vev to cancel the bare p vac and give 
P^ c ~(10- 3 eV) 4 . 

The model has four adjustable couplings (k, J, £ and p vac ) and accounts 
for two observed quantities (k b s and H). So we can express two of our 
couplings - say k and J - in terms of £, p vac and the measured parameters 
K eS , if and pH = 3(mf) 2 H 2 . 

(mp) 2 = 1/k 2 = (p vac + 3p£) (24) 

J 2 = 12^ 2 (p vac - pfj (25) 
The scalar vev can likewise be expressed as 

02 = 12^ ( Pvac ~ P ™ c) • (26) 

Note that (25) requires £p vac > Just to plug in some numbers, if we take 
£ = 1 and p vac = (10 18 GeV) 4 we would have 

J ~ (10 MeV) 3 m P « (0) ~ 10 86 eV . (27) 



7 Although we are most interested in having £ and p vac positive, there are consistent 
solutions in which £, p vac and the bare k 2 are all negative. 
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A few comments on these results: 



1. The see-saw mechanism requires J 2 <C £ft 2 p 2 ac . Here however we are 
interested in the regime p vac ^> p^f c . In this case (24) and (25) imply 



J ~ C K Pvac; so m ^ ac ^ the see-saw mechanism isn't operative. Rather 
we have an approximate cancellation of the vacuum energy, p°f c <C p vac , 
but also an approximate cancellation of the Planck mass, trip -C mp, 



such that curvatures of the two solutions in (|20|) are roughly the same. 

j2 ^ t„1 2 

vac 



Indeed for J 2 ps £k 2 p 2 (20) reduces to 



1 2 ^ . 1 



The see-saw regime of the i?0 2 model is actually rather unphysical, as 
one can show that the low-curvature see-saw solution has both p^ c < 
and (mf ) 2 < 0fl 

2. The model can be made compatible with solar system tests of gravity. 
This is unlike the original R+l/ R models of cosmology which, if naively 
extended to the solar system, ha ve be en ruled out [HI [TTJ . As we 



discuss in more detail in appendix A.2 the R<p model can be related 



to Brans-Dicke theory, and solar tests require 



eff 

S I Pvac 



vac/ 



With a natural value for p vac this means £ must be extremely small, 
perhaps |£| < 1CT 125 ; the sign of £ is correlated with the sign of p vac . 
Then the vev of <fi becomes much larger than the Planck mass, ~ 
mp/v^, which raises an important issue: our conclusions are sensitive 
to the presence of additional non-renormalizable operators suppressed 
by powers of (p/m P . 



3. As we show in appendix A.l the low-curvature de Sitter solution can 
be made metastable over cosmological timescales. Interestingly, the re- 
sulting requirement on £ is similar to that from solar system constraints. 
Therefore, according to this model, if the dark energy equation of state 
is observed to differ from —1, one should also expect to see general 
relativity violations just around the corner in solar system tests. 



5 We are grateful to Ali Masoumi for pointing out this pathology. 
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4. Achieving p^f c -C p vac and p v f c <C (mp ) 4 requires severe fine-tuning. 



In particular (24) can be rewritten as 



Pvac ^ _ Pvac . Pvac _ ^-240 (2g) 



where we've assumed p vac ~ {m e p) A . Likewise (24) and (25) can be 
combined to give 

J ^pL -1 + 2^ = 1 + O(10- 120 ) . (30) 

S Pvac Pvac 

These two conditions mean the Rcf 2 model is actually much more fine- 
tuned than a pure cosmological constant. This is in accord with our 
general discussion of tuning in F(R) models at the end of section |3j 



As we show in appendix A. 4 the necessary tunings are destabilized by 
radiative corrections. The easiest way to see this is to map the model 
to Einstein frame, where no symmetry protects the vacuum energy. 

rrfp sets the scale at which gravity becomes strongly coupled, so one 
might argue that the natural scale for vacuum energy is (m c p) 4 . This is 



the value we adopted in (27). It leads to an amusing coincidence: with 
Pvac ~ ( m< p) the total energy inside our Hubble volume, estimated 
as pl^ c /H 3 , is always comparable to the bare Planck mass nip. One 
could very well criticize this choice, on the grounds that p V ac can't 
know about ^-dependent quantities such as mfF; this is related to the 
question of why the potential for is so flat. To address this one could 
imagine tying p vac to a non-gravitational scale such as the scale for 
supersymmetry breaking. 



5 Gauss-Bonnet model 



The mechanism we have discussed is clearly rather general: one could start 
with a model in which vacuum energy is redundant and explicitly break 
the redundancy in any number of ways. One possibility is to introduce a 
scalar field with a non-minimal coupling to a more general curvature invari- 
ant. Upon integrating out the scalar, such a model would yield one of the 
generalized modified gravity theories studied in [IS] . 
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A particularly appealing possibility is to couple the scalar field to the 
Euler density Q . Thus we consider the action 

S = J d A x^~g (J^R - l -d^d^ - l -in 2 G<p 2 + J0 - p vac ) (31) 

where 

g = R 2 - AR^ v + R^xaR^ ■ 

This is known as Gauss-Bonnet gravity. For a recent review see [19J. Gauss- 
Bonnet gravity has a number of desirable features. It gives a ghost-free 
theory with second-order equations of motion (20]. It can be generated by 
a 1 corrections in string theory |21j. And most importantly from our point of 
view, it will allow us to avoid many of the fine-tuning problems of the Rip 2 
model. One could imagine adding a mass term for the scalar, but the mass 
would have to be tiny: as can be seen below, we need m 2 < £,k 2 H 4 for the 
model to work. 

Assuming is constant we can integrate it out, <fi = J/^k 2 Q to obtain 
the effective action 

This is a particular example of the generalized modified gravity of Carroll 
et. al. [TS], supplemented with a vacuum energy term. On maximally- 
symmetric spaces Q = R 2 /6. So not surprisingly there are solutions in 
which the Gauss-Bonnet term dominates over Einstein-Hilbert and R 2 ~ 

J 2 / £,K 2 Pvac- 



To analyze the model more carefully we should vary (31) with respect 
to cf) and the metric then restrict to constant <ft an d maximally symmetric 
metrics. We do this in appendix [Cj However if one is only interested in de 
Sitter solutions a general analysis isn't necessary. Recall that Q is a total 
derivative in four dimensions, so upon plugging a constant vev for <p into 



the action (31) the Gauss-Bonnet term can be dropped. The effective Planck 
mass is therefore unshifted from its bare value, nfp — trip, while the effective 
vacuum energy is pl^ c = p vac — J(<f>). Thus the Friedmann equation is 

R = 12H 2 = 4K 2 pll = An 2 (p vac - J<P) (32) 

where the scalar equation of motion fixes 

= QJ/£ K 2 R 2 . (33) 
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Eliminating gives a cubic equation for R. The see-saw mechanism operates 
when J 2 <C £ K Vvac- To leading order in this regime the three solutions are 



a 72 

i^^pvac, R^±J— . (34) 



vac 



We require that £ and p vac have the same sign so that all three roots are real. 
Then we have one solution with large curvature and two solutions with small 
opposite-sign curvatures. 

In the small-curvature de Sitter solution we can solve for J in terms of £, 
p vac and the observed quantities k and H. 

J 2 = 24£ K 2 # 4 pvac (35) 

Likewise the vev is given in terms of these quantities by 

(0 2 > =p vac /24^ 2 # 4 . (36) 

Just to plug in some numbers, setting £ = 1 and p vac = (10 18 GeV) 4 we have 

J~(l(T 13 eV) 3 (0)~lO 146 eV. 

A few comments on this model: 

1. Unlike the Rcf 2 model, the Gauss-Bonnet model realizes the see-saw 
mechanism in the regime of interest: an underlying Planck-scale vac- 
uum energy can drive either a fast or a slow de Sitter expansion. 

2. There is no need to fine-tune the bare Planck mass, nor is there any 
need to take the bare p vac <C rrip. In this sense the Gauss-Bonnet model 
avoids some of the fine-tunings necessary in the Rcj) 2 model. 

3. The only small number we need in Gauss-Bonnet is J 2 -C £ft 6 p 3 ac . Note 



that (|35|) implies 

T 2 m 6 / n eS \ 2 

p - ' Pvac » - itr 240 . (37) 



£Pvac VP 



As we will show in section 5.2 this small parameter is stable under 
radiative corrections! 

The Gauss-Bonnet model does, however, involve super-Planckian vevs. 
This makes the model sensitive to any other higher-dimension operators 
that might be present in the Lagrangian. 
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5.1 Classical stability and solar system tests 



In appendix [C] we show that the low-curvature de Sitter solution is metastable 
over cosmological timescales provided £ < 1. However we also want the model 
to be compatible with solar system tests. Without performing a full analysis, 
we believe that for sufficiently small £ the model will pass all solar system 
tests of gravity. The essential point is that the scalar field can be made 



arbitrarily weakly coupled. To see this take the Gauss-Bonnet action (31) 
and set 

(f) = ((f)) + 5(f) 9ilv = ( g ^) + Sg^ . (38) 

Dropping a total derivative, the full Gauss-Bonnet action becomes (here p°f c 
is a fixed numerical quantity, and R for example means (R) + 5R) 

S = J d^x^f-^R-p^-^d^d^)-^ 2 ^) <^-^ 2 S(^) 2 

(39) 

The last two terms represent non-minimal couplings between 5<f> and the 



metric. But for fixed p vac , (36) implies that ((f)) ~ So as £ — > these 

non-minimal couplings vanish (they scale as £ 1//2 and £, respectively). Thus 
as £ — > the model goes over to Einstein gravity with a cosmological constant 
plus a massless, minimally-coupled scalar. Such a model is compatible with 
all solar system tests of gravity, so we expect Gauss-Bonnet to pass solar 
tests provided £ is sufficiently small. This does not address the question of 
exactly how small £ must be; to settle this question would require a more 
detailed analysis along the lines of [22J.H 



5.2 Radiative corrections 



Finally we consider radiative corrections in the Gauss-Bonnet model. We 
will argue that radiative corrections to the scalar potential are under control 
provided £ is sufficiently small, roughly £ < 10~ 240 . In this sense, for small £, 
the Gauss-Bonnet model provides a technically natural explanation for dark 
energy. 

9 The analysis in [25] has a restricted range of validity, see their (20), and cannot be 
directly applied to our model. 
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When £ = J = the scalar field has a shift symmetry 0-^0 + const, 
which forbids any corrections to the scalar potential. So one might expect 
radiative corrections to the scalar mass, for example, to be proportional to 
£. On dimensional grounds one might expect a correction 5m 2 ~ t° 
be generated. If the UV cutoff scale Ajjv ~ this would be bad news, 
because the Gauss-Bonnet model requires m 2 < C,k, 2 H 4 . 

Here the remarkable structure of Gauss-Bonnet gravity comes to the res- 
cue. Expanding about a de Sitter solution, the relevant interaction vertices 



can be read off from the last two terms in (39). Since the Euler density is 



a total derivative, the vertices vanish when all scalar lines together carry no 



net momentum into any vertex 10 This means graphs like (scalar lines are 
solid, graviton lines are dotted) 



which naively generate an correction to the scalar mass, actually vanish 
at zero external momentum. The leading correction to the scalar mass seems 
to come from a two-loop diagram. 



A simple estimate is that this diagram generates an 0(^ 2 ) correctiorj"] 
x 2 <-2 4 f d 4 p d 4 q 4 4 k 2 k 2 1 



(2tt) 4 (2tt) 4 * * p 2 q 2 (p + q) 2 

e^ v (4i) 



10 Additional interaction vertices arise from the scalar kinetic term, but due to the shift 
symmetry they vanish at zero scalar momentum and we can ignore them. 

n Each vertex has a coupling £k 2 and involves two curvatures, hence two powers of 
momentum for each graviton. Setting g^ v — (g^ + Sg^ the graviton propagator ~ n /p . 
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With Auv ~ m P this means Sm 2 ~ £, 2 m P . So radiative corrections to the 
scalar mass are under control for 

iVp < ^ 2 H 4 

or equivalently 

£<(pf c /- 4 p) 2 ~ 10- 240 . 
In a similar way one can study corrections to the linear source J. Here 



preserving the relation (37) requires that 5 J/ J < 1. As above, the leading 



correction to J seems to come from a two-loop tadpole. 



This generates a correction 
5J 



2 . > / >, / d^p d 4 q 4 a k 2 k 2 1 

& m 



v q 



(2vr) 4 (2vr) 4 ^ p 2 q 2 (p + q) 2 



£V(0>A 10 



uv 



(42) 
(43) 



Taking A;yy ~ rap, and using (33) to eliminate the vev, this means 



5J C m p 
J H A ' 

Again radiative corrections are under control for £ < 10~ 240 . 



Finally, what about radiative corrections to p^ 
tive corrections to generate a shift 



itself? We expect radia- 



A 4 



■> (rap) 4 , this means 5p vac ~ p vac - This preserves 
the tuning (37), so the model is stable under radiative corrections to the 



Assuming we start with p v 



vacuum energy. More directly, as can be seen from (32) and (34), the model 
will compensate for the change in p vac by shifting (0) so as to preserve a small 
effective vacuum energy, with Sp^ c ~ p^f c . For this argument to hold it's 
important that shifting (0) doesn't change the Planck mass. (In this respect 
the Gauss-Bonnet model differs from the R<fr 2 model studied in appendix 
[Ol) 
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6 Conclusions 



In this paper we explored the idea that the underlying vacuum energy could 
be large, of order the (effective) Planck scale, with the observed slow Hubble 
expansion due to the dynamics of a scalar field with non-minimal couplings. 
We showed that the idea could be realized in the R<p 2 model and its F(R) 
generalizations. However all such theories, involving only the scalar curva- 
ture, require severe fine-tuning and are unstable with respect to radiative 
corrections. 

We went on to consider the Gauss-Bonnet model in which the scalar field 
couples to the Euler density. The Gauss-Bonnet model requires that one 
combination of parameters be small, namely J 2 <C £>*Vvac- However we 
argued that for small £ this parameter is stable under radiative corrections, 
and in this sense the Gauss-Bonnet model provides a technically natural 
explanation for dark energy. 

The crucial feature that makes all this possible is the fact that the bare 
vacuum energy is a redundant coupling, and therefore unobservable, when 
£ = 0. This allows us to shuffle the necessary fine-tunings among various 
couplings in the Lagrangian. The way in which tunings are shuffled will 
play an important role in determining how likely these models are in the 
landscape. Also, by shifting some of the tuning to the Gauss-Bonnet term, 
pseudo-redundancy made it possible to construct radiatively stable models 
for dark energy. 

This leaves many open questions: 

• Can one envisage a cosmological history in which the scalar field nat- 
urally evolves to take on its required expectation value? Or does a 
realistic cosmology require modifying the model in some way? If one 
does keep £ small throughout cosmological evolution then the quantum 
diffusion studied in [231 121] will play an important role. 

• Some models, in particular Gauss-Bonnet, realize a see-saw mecha- 
nism in which a large underlying vacuum energy can drive either a fast 
or a slow de Sitter expansion. Can one find interesting cosmological 
solutions in which the universe spends some time inflating near the 
high-curvature solution before evolving to low curvature? 
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• The models we studied require super-Planckian vevs which are often 
viewed as problematic [25]. But are these large vevs really necessary? 
Or could the models be modified in some way to eliminate them? 

• Can the necessary non-minimal couplings be realized in a UV-complete 
theory such as string theory? Or is there some fundamental obstacle 
to achieving this? 

• One intriguing feature of the Rcf 2 model is that cosmological stability 
and solar system constraints put roughly similar requirements on the 
coupling £. Therefore, if £ takes a value such that the dark energy 
equation of state deviates from —1 by an observable amount, it is likely 
that one would also see violations of general relativity in precision solar 
system measurements. This offers an interesting way to test the model. 
It would be useful to check if the same property holds in the Gauss- 
Bonnet model. 

Acknowledgements 

We are grateful to Matt Kleban, Ali Masoumi, Alberto Nicolis, Massimo 
Porrati and Iggy Sawicki for valuable discussions. PB, KH, LH and DK are 
supported by DOE grant DE-FG02-92ER40699 and by a Columbia Univer- 
sity Initiatives in Science and Engineering grant. 



A F(R) dynamics 

In this appendix we study the dynamics of the F(R) type models in more 
detail. We study classical stability of the de Sitter solutions, discuss solar 
system and cosmological constraints, and analyze radiative corrections, in 
some cases in general F(R) gravity, and in some cases in the context of the 
R<j) 2 model. 
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A.l Classical stability 



To study the stability of the F(R) solutions it's useful to make a conformal 
transformation. For a generic scalar-tensor theory in Jordan frame 

S = J d 4 x^9 - \h($)g» v d^dA - 

we define = 2(k c r) 2 f(4>)g^u, where K e ff is a fixed numerical quantity, to 
go to the Einstein-frame action 

S = J d A x y^g | Imr^d^-ViA . (44) 

In this frame the scalar field is minimally coupled and 

2{k cS ) 2 \f f 2 j (2(K eS ) 2 f) 

By the definition of the effective Planck mass note that = l/2(K e s) 2 . 
Thus on our static solution the Jordan-frame and Einstein-frame metrics are 
the same. We can then read off the dynamics of the scalar field just from 
its kinetic term and potential. The key point is that by making h\u\ large 
we can suppress any temporal variation of the field and make the solutions 
(meta-) stable over an arbitrarily long time scale. 



For further discussion let's specialize to R<p 2 model for which^ 

f{4>) = _ tftf) h {cj)) = i v(<t>) = Pvac - J0 . 

In the regime of interest, where |£| <C 1, we also have h\u\ ~ 1. So small 
fluctuations about the static solution are governed by the effective potential 



f k Y p 



V= —) 



vac 



2A2\2 



The potential is sketched in Fig. 1. It vanishes when = p vac /J and blows 
up where the effective gravitational coupling diverges, at <fi 2 = 1/^k 2 . The 
potential has two critical points, with large and small vevs, corresponding to 
the two solutions in (20). We are interested in the large- vev solution. 



12 By a field redefinition, instead of making h large, we will set h — 1 and take £ to be 
small. 
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V(<f>) 




Figure 1: Effective potential for the R<fr 2 model. 



To study stability we expand about the unstable point, (f> = {(f)) + 5(f). To 
quadratic order the potential is (for p^ c <C p vac ) 



r vac 



^Pv 



[m 



cS\2 



+ omf 



The timescale associated with the instability is 



—V" 



'1/2 



eff xl/2 
rvac \ 



6^Pv 



So for sufficiently small £, say £ < 1CT 120 , the solution will be stable for more 
than a Hubble time (and Hubble friction will play an important role in the 
evolution of 6). 



A. 2 Solar system tests 



In section [3] we constructed F(R) models that give the right Planck length 
and de Sitter radius. However we also want models that are compatible with 
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solar system tests of gravity. One way to achieve this is to give the scalar 
field a non- minimal kinetic term. As in (|8|) we take 



S 



d 4 x^(f(<p)R-h(<p)d, 



- v{<p) + s_ 



We assume that matter is minimally coupled to g^ v . By taking h large any 
spatial variation in is suppressed and the model can be made compatible 
with solar system tests of gravity. Note that by a field redefinition <fi — > 
(j)/ a/ (h) taking h large is equivalent to making the scalar field very weakly 
coupled, i.e. £ small in the Rep 2 context. 

To show this in more detail we use the equivalence to Brans-Dicke gravity 
[26J. The classic solar system constraints on Brans-Dicke gravity were ob- 
tained by setting V = 0. We can therefore translate the known constraints 
to our situation as long as terms related to V in both the gravitational and 
scalar equations of motion are subdominant in the solar system. The Einstein 
equations read 



G m = r 1 {hT^ + ±1% + V M V„/ - g^Uf 



(46) 



and the scalar equation of motion is 



, ^ , ldh„ u , dv d f r, 



(47) 



^From these equations it seems we need 

df 



fR^V , 



dV 



(48) 



where R is the curvature in the solar system. Recall that our static de Sitter 
solution has foRo = 2Vo and f^Ro = Vq, where R ~ Hq -C -R so iar, and f and 
Vq denote cosmological vevs. It therefore seems the above two inequalities are 
easily satisfied. However there is an implicit assumption here, that the value 
of <j) in the solar system does not differ too significantly from its cosmological 
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value, so that / so iar ~ fo and V^ i ar ~ Vo- This will be true for sufficiently 
large h, as we now showp*] 

To study the variation of (f> within the solar system we linearize the scalar 
equation of motion, setting cf> = ((f)) + 5(f) with 5(f> obeying 

V 2 5<p + f ° R - V U <p + m ~ Ro) = . (49) 

For large ho we can neglect the mass term. The source term is dominated by 
the scalar curvature in the solar interior, where R ~ ^ ff p . Thus to a good 
approximation 

V 2 ^ + / ° K f M0 ^(x)=O. (50) 
So 5(f) is maximized at the surface of the sun, where 

47rh r Q 

^From @ and Q it follows that 5V/V « 2 5///, so we only need to impose 

y ~ 4nhoV r Q K< 1 " (52) 

This is easily satisfied for large ho- Then the model is equivalent to Brans- 
Dicke theory, and solar system tests mostly bound two parameters j27] 



7= = l±O(10" 5 ) , 5=1 + 7 ^77 r = 1 + 0(10 

' 2 + cu v ; ' M (3 + 2cj) 2 (4 + 2w) v 

where uo and u/ are related to our / and h by: 



(53) 



"-HIT (54) 

' = 2 f ( ldh f + 1 h _ r ^ 

7 V2#(d//#) 3 2{df/d<f)) 2 {df/d<t>) A )' { ' 



13 As an alternative approach to satisfying solar tests, imagine taking /i = 0. Then the 
model is strictly equivalent to F(R) gravity, even within the solar system. The scalar 
equation of motion R = V / f requires that the scalar field, or at least V /f, change 
significantly inside the solar system. Then with suitable F(R) one can evade solar tests 
of gravity via an analog of the chameleon mechanism [15j. 
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As long as h is large enough the solar system tests are passed. 



For example, as shown in [26J, in the R<p 2 model we have h = 1, / = 
- I^ 2 , V = p vac - J<p and therefore u = With the help of 



2^2 - V = Pvac - 

this can be rewritten as 



off 
r vac 



Hp 



vac Hva,c) 



(56) 



The bound (53) requires uj > 10 5 (we take uo > —3/2 so that in Einstein 



frame the Brans-Dicke scalar has a right-sign kinetic term). This leads to 
the constraint on £ in ( 28 ) . However for this constraint to be valid we should 
make sure that ( 52 ) is satisfied. Using J(f> w p vac ( 52 ) can be rewritten as 



cflf 
r vac 



^0 



4tt • 10" 120 • 10 5 



(57) 



This is less stringent than (28) 



A. 3 Matter and radiation domination 

The de Sitter solutions we found give us phenomenologically acceptable cos- 
mic acceleration today. Here we study the conditions under which these 
solutions can be extended to the early universe in the presence of matter and 
radiation. More concretely we seek a solution of the form = 0o + <f>i where 
4>o is a constant chosen to give us phenomenologically acceptable cosmic ac- 
celeration today and 0i is a small time-dependent perturbation. It is by no 
means obvious that this is the most interesting dynamical solution, but such 
a quasi-static solution for the scalar field, if it exists, provides a particularly 
simple way to satisfy all observational constraints. 

Applying the FRW ansatz to the equations of motion of the Rcf 2 model 
we find 



-R 



[ \(-P + 3p) mattcr + \{H 2 - 4(p vac + J0)) (58) 
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where we have taken the trace of the Einstein equations (46). The scalar 
equation of motion is 



h{4> + 3H4>) + -h'<p 2 = - J - £(f>R . 



(59) 



The static de Sitter solution obeys 

J + &0R0 = , Ro = 4(p vac + J<po)/(M 2 - Z<Pl) 
where Rq = 12ifo is the curvature today and where we defined 

Pvac = Pvac + J0, Mj = M* - ^ . 



(60) 



(61) 



Here is the bare Planck mass (the quantity we denoted mp previously). 

Going back in time, the source term for the (f) equation is J + £,<pR ~ 
£,4>R because R oc a~ 3 during matter domination (we will discuss radiation 
domination below). Assuming (f>i <C <p we therefore have 



H d 

a 2 da 



da 



n 



(62) 



where we have used 0i + 3H<px = (H/a 2 )(d/da)[a A Hd<j)/da\ and R = 3H 2 
during matter domination. Integrating we have 



f da 






J a 2 Hda 




1 a A H 





da 



+ B 



(63) 



where A and B are dictated by initial conditions. The constant B can be 
absorbed into 0q. The mode associated with A decays with time, and one 
can reasonably argue that for a wide range of initial conditions this mode 
will be subdominant at all times relevant for observation. There remains the 
driven mode, given by 



-2—00 l n cl . 
h 



(64) 



Clearly one can make <p\ <C <po by choosing £/h <C 1. The weak dependence 
on a means it is easy to satisfy <pi <C <pQ for all observationally relevant 
scale factors. However to be self-consistent we also need to check that the 
perturbation does not dominate over the matter energy density and pressure 



25 



in (58). Examining the various contributions from (f>i and its derivatives, 
the dominant contributions are of order (£ 2 /h)H 2 <pQ (such as from tfi 2 ), and 
the self-consistency requirement is that this is much smaller than 

pmatter ^ 

H 2 M 2 V In other words, we need to make sure that (6//i)£</>o < M pi- Recall 
that the phenomenologically interesting static de Sitter solution has M 2 3> 



M p 2 ,, 



which basically means £0q ~ M 2 (see (61)). Therefore a solution for 
the scalar field that does not significantly moony the expansion rate at early 
times requires 



S „ pi rvac 

— <<C^ rs_/ 

h M 2 p vac 



(65) 



This requirement is rather severe, but it happens to be about the same as 
the requirement for stability of the static solution. 



What happens if we go back even further, to the era of radiation domi- 
nation? An interesting observation is that due to its traceless stress tensor 
radiation makes no contribution to R. However this doesn't mean R van- 
ishes, due to the stress tensor for 0, so in fact a quasi-static scalar field would 
mean R ~ Rq + SR during radiation domination. Then the <fi equation of 



motion (59) tells us 



n 



h 

R (, 



,5R 



hM 2 , 



R 



rh W° 

h 1M 2 



XMl 



hM 2 pl 



pi 



KMl 



+ 3#0i) 



3#0i) 



(66) 



where we've used the fact that J + £0o-Ro 



0. 



We're allowing 0i to evolve on cosmological time-scales, so H<pi ~ H 2 (pi 
which is much larger than (^//i)(M fe 2 /Mp 1 )i?o0i 14 This means we can drop 



the first term on the right hand side of (66). ATso given (65) we can drop 



the last term on the right hand side of (66). So we're left with a very simple 
equation 



+ 3#0i ~ 



(67) 



14t 



From (651 we have (£//i)(Af fc 2 /M^Rofa < i? </>i 



H 
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with solution 



(p l = A/a + B. (68) 

A and B are determined by initial conditions. B can be absorbed into 0o> 
and A can always be chosen to be small enough to satisfy any self-consistency 
requirements one needs to impose. (Think of it this way: suppose A/a takes 
some natural value in the very early universe, then for the a's of interest it 
is likely that A/ a is very small.) The upshot is that a quasi-static approxi- 
mation for the scalar field is valid during radiation domination as well. 

In summary, the presence of matter and radiation is consistent with a 



quasi-static scalar field as long as (65) is satisfied. A similar analysis with 



similar conclusions can be made for the Gauss-Bonnet model. 



A. 4 Radiative stability 

So far we have treated the F(R) models as classical field theories. But one 
of the central mysteries of the cosmological constant is whether it can be 
protected from large radiative corrections. In this section we study quantum 
effects. To be concrete we will focus on the R<p 2 model, although similar 
results should hold in general F(R) models. 

First let's consider the model in Jordan frame. 

A key observation is that £ = J = is an enhanced symmetry point, where 
the theory is invariant under shifts — > <fi + const. This means radiative 
corrections to quantities like £, J and the scalar mass m (set to zero above) 
have to be proportional to the symmetry-breaking parameters £ and J them- 
selves. To summarize the results of our analysis: radiative corrections to 
these symmetry-breaking parameters are under control if £ is sufficiently 
small, roughly |£| < 1CT 120 . However this argument does not address radia- 
tive corrections to the vacuum energy p vac itself. Indeed we will find that, 
from the point of view of quantum corrections to p vac , the Rcj) 2 model is just 
as unstable as a bare cosmological constant. 
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We will perform the rest of our analysis in Einstein frame. 

s = Jd A xy^ - Imr^d^ - v(4>) 

We will concentrate on quantum corrections to the scalar potential V . Before 
studying our model, consider a scalar field theory in flat space with a potential 

= P + J<t> + \rn 2 <p 2 + + ^A0 4 . 
In this flat-space theory a one-loop vacuum bubble should make a correction 

where kuv is a UV cutoff scale. A one-loop wart diagram should make a 
correction to the scalar mass 

=>- Am 2 ~ XAfjy . 






Finally a one-loop tadpole diagram should make a correction to the source 

AJ ~ gA 2 uv . 

This will generate a shift in the vev 

A(<f>) = -AJ/m 2 
which in turn produces a shift in the vacuum energy (for J = 0) 

A(V) = (AJ) 2 /2m 2 . 

What do these results mean for the R(p 2 model? Setting = (0) + 5(f), the 
effective potential of the R<p 2 model has an expansion in powers of J5<p/ p^Ff] 



15 Here p°f c is a fixed quantity, independent of <fr. It's useful to note that 

J 2 _ 4g(Keff) 2 j0vac 

(„cff \2 - „off 
Vrvac/ Fvac 
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There are three classes of corrections we will study. 
Corrections to m 2 . 



In Einstein frame the scalar mass i; 



16 



^Einstein ~ J / Pvac • 

Given the quartic coupling in our model the radiative correction gives 

A™lin S tein ~ J*tijvHj&? 

or 

^ m Einstein J ^UV ^Pvac A uy 



m 2 (rt eS ) 2 n eS (m cS ) 2 ' 

'""Einstein iPvac/ Fvac \" L P I 

Radiative corrections to the scalar mass are stabilized for £ < 10~ 120 , assum- 
ing p vac ~ (infY and A uv ~ mf. 

Linear source term 

In Einstein frame we'll generate a linear source term for the scalar field 

A J ~ J 3 A-uv/ (P^L) 2 ■ 
This shifts the vev and in turn corrects the vacuum energy. 

AV ~ (AJ) 2 /m| instein ~ J*Alv/(f£f 

The quantity we need to look at is 



AV J 4 A UV /£PvacV A 4 



Pf c (Pfj 4 \P e JLJ {mf* 



uv 



vac Vrvac/ \ rvac / v'^P 

Again radiative corrections are stabilized for £ < 1CT 120 , assuming p vac ~ 
{m e p) A and Auv ~ mfF. 

16 In Jordan frame the scalar mass is much smaller, mj ordan ~ ^H 2 ~ C( K eff) 2 Pvac ~ 
J 2 / Pv&c- This fits with our stability analysis: given the Jordan-frame mass you might have 
thought £ = 0(1) would be sufficient to guarantee stability on cosmological timescalcs, 
but in Einstein frame you find that you actually need £ < 1CP 120 . 
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Corrections to p' 



.eff 

vac 



Here the R<p 2 model suffers from the usual instability of vacuum energy 
with respect to radiative corrections: in Einstein frame, no symmetry forbids 
graviton loops or loops from generating 

APvac ~ Kv ■ 

This means the R(p 2 model doesn't avoid the usual naturalness problems 
associated with vacuum energy. Although this is easiest to see in Einstein 
frame, one can reach the same conclusion in Jordan frame. In Jordan frame 
graviton and loops, or loops of standard model particle fields, should gener- 
ate a shift Ap vac ~ h^jy. This change in p V ac can be compensated by shifting 
(4>) to preserve a small p v f c , but the necessary change in ((f)) completely desta- 
bilizes the effective Planck massj^] To be more explicit about the difficulty, 
suppose we have a correction in Jordan frame 

Ap vac ~ Afjy . 

You might think we've achieved something: with p vac ~ (m e p) A and Ajjy ~ 
m c p, don't we naturally have Ap vac /p vac = (9(1)? Unfortunately the relation 



(30) 



^4 = i + io- 120 

>r vac 

requires Ap vac /p vac < 10~ 120 . Equivalently, it requires Ap vac < p v f c . So the 
Rcf) 2 model is just as sensitive to radiative corrections as a pure cosmological 
constant. 

To summarize, with £ < 10 -120 radiative corrections to the non-constant 
parts of the scalar potential seem to be under control. It's curious that 
this value of £ was already required for cosmological stability. However the 
constant term in the potential is not protected, and makes the R<p 2 model 
just as unstable with respect to quantum corrections as a bare cosmological 
constant. We expect this last conclusion to apply to any F(R) model, since 
any F(R) model can be given a scalar-tensor description and mapped to 



Einstein frame as in (44). 



17 The Gauss-Bonnet model gets around this because a change in (0) does not affect the 
Planck mass. 
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B Further examples of F(R) 



In this appendix we present further examples of modified gravity actions, 
built purely from the scalar curvature, which have low-curvature de Sitter 
solutions. As in section [3] we will write the action as 

S = J d*xV=g(mf fF(y) 

where 



y = R/A{mf 



eff\2 



Recall that any function F(y) satisfying (16), namely 



F{e) = e F'(e) = 2, (70) 

will have a low-curvature de Sitter solution with the right effective Planck 
mass, where e = 3Hl/{m c p) 2 ~ 1CT 120 . 



Clearly an infinite number of functions satisfy ( 70 ) : F(y) = exp (y / e) — 1 



1 + 2y 2 /e — 1, 2e \og(y/e) + e, ... To be a bit more systematic about this 
we make the ansatz 

oo 

F(y) = a -y n ■ ( 71 ) 



n=— oo 



Then ( 70 ) implies the conditions 



Here are a few examples. 

Example 1. We can satisfy the necessary conditions with a single cou- 
pling constant by taking F(y) = y 2 /e. That is, the gravity action 

S = [ d A x^R 2 /lQe (73) 



has a solution with both the correct Planck mass and the correct vacuum 
energy. To make the model compatible with solar tests we invert the Legendre 
transform and add a large scalar kinetic term. 

S = J d 4 x (Vr - h(d4>) 2 - 4e0 2 ^) (74) 
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One curious feature: neglecting the scalar kinetic term, this action is invariant 
under the scale transformation 



g^u -> ag^ <p—><p/a. (75) 

This means the model has de Sitter solutions with arbitrary radius (and 
varying Planck mass). It is an extreme example of the fact that, although 



(70) guarantees the existence of a solution with the desired properties, it does 



not preclude other de Sitter solutions. 

Example 2. With two coupling constants, the most obvious way to 



satisfy (70) is to take F{y) = 2y — e. This just corresponds to ordinary 
Einstein gravity with a small cosmological constant. Another option is to 
choose F = —2e 2 /y + 3e. Going to the scalar-tensor description, and adding 
a large scalar kinetic term, the latter example shows that an Einstein-Hilbert 
term is not necessary to pass solar system tests. 

The examples we have considered so far do not interest us much, because we 
want models with do ~ — 1 representing a Planck-scale bare vacuum energy. 
Our remaining examples have arbitrary a$. 

Example 3. Choose 

a = arbitrary , a\ = — - — - , a_i = — e(a + e)/2 (76) 

and the rest of the a n 's vanish. That is, the gravitational action is 

S = (,nfy J d'x^g + + a ) . (77) 

Inverting the Legendre transformation, and choosing the relation between / 
and 4> to set (p = (m e p) 3 /R, this is equivalent to the R<p 2 model studied in 
section HI The scalar-tensor action is 



S = d 4 X\f^g 



ai(m c p) 2 



Aa^R-hid^) 2 -V((j>) 



V(4>) = -8(mf ) 3 a_!0 - {mfYa (78) 
where we added a kinetic term for (p. 

Example 4. Choose 
a = arbitrary , a_i = — 2e(a — 2e) , a_2 = e 2 (ao — 3e) (79) 
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and the rest of the a n 's vanish. That is, the gravitational action is 



«) 4 



4(mf) 2 a_ 1 16(mf) 4 a_ 2 



R 



+ 



R 2 



+ a 



(80) 



This looks strange, but remember that for solar system tests one must revert 
back to the scalar-tensor theory. Inverting the Legendre transform, and again 
taking = (m e p) 3 /R, we have the equivalent scalar-tensor theory 



V{<j>) = -8(m^) 3 a_!0 - 48(m^?) 2 a_ 2 2 - (mf) 4 a 



(81) 



C More on Gauss-Bonnet 



In this appendix we study the classical stability of solutions to the Gauss- 
Bonnet equations of motion, and show that for £ < 1 the low-curvature de 
Sitter solution is metastable for more than a Hubble time. 



C.l Equations of motion 



We consider a generalized Gauss-Bonnet action coupled to matter, 



s 



M 



where the Euler density is 

G = R 2 - AR^R, U + R^ kX R^ kX . 
Define the following tensor 
1 5 



-9&9. 



= 2RV»V v f{4>) + 2g^ [2^V^V7(^) - RDf ((/>)] 
-8 J R^V" ) V p /(0) + 4/TTJ/(0) - 4 J R ww V /0 V CT /(</>). 



(82) 



(83) 



(84) 
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The equations of motion for the metric are then 



1- 



2 

where 



7^ + 7^-2^ 



(85) 



1 



T+, = h ^ ~ 2^(^) 2 ) " 9^V{<j>) . (86) 

The equation of motion for the scalar field is 

/i(0)D0 + h'^)d^d,4> - v'(4>) - f(4>)g = o . (87) 

It's a remarkable property of Gauss-Bonnet gravity that one obtains second- 
order equations of motion, even though one is starting from a higher-derivative 
action which naively would be expected to generate fourth-order equations 

In maximally symmetric spacetimes we have 



Rfiu ^RQfiu ^RQijlv 

q = r 2 /6 

0^ = 0. 



In the absence of matter, and assuming a maximally symmetric metric and 
constant 0, the equations of motion simplify to 

R = 4kV(</>) 

v^) = -\f^)ie. (88) 



These are the Gauss-Bonnet analogs of (J6|, Q in F(R) gravity. We see that 
critical points of the scalar potential are not what give us de Sitter solutions. 
Rather we must look for points on the potential such that 

H0) = -^ 4 f(0M0) 2 . (89) 
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C.2 Classical stability 



Here we show that for £ < 1 the low-curvature de Sitter solution is metastable 
on cosmological timescales. Plugging in the FRW ansatz 



ds 2 



-dt 2 + a 2 (t)\dx\ 2 , = 0(t) 



(90) 



we obtain the Friedmann equation from the pv = 00 component of (85) 

h 



K 



-0 2 - V{4>) - 240f (0)# 3 = pu 



(91) 



where pu is the matter energy density. Assuming constant h, the scalar 



equation of motion (87) becomes 



h (0 + 3#0) + F'(0) + 24f(0)# 2 [H + H 2 ) = . 



(92) 



In these equations H = a/a is the Hubble parameter and some useful iden- 
tities are 

R = QH + 12H 2 , g = 2AH 2 {H + H 2 ) . (93) 



rather are consequences of (91) and (92) 



The remaining components of (85) do not yield independent equations, but 



We are interested in fluctuations about a low curvature de Sitter solution 
in which H and are constant. Setting 



H(t) = H + SH(t), 0(t) = 0o + 6<f>(t) 



to first order in the fluctuations (91) reads (in the absence of matter) 



6H 

24/'# o 3 <50 + V'5(j) ^5H = 



(94) 



(95) 



where /' and V are evaluated on O . Similarly (92) reads 



h5(j> + 3hH 5<j) + (V" + 24f"H%) 50 + 24f'H 2 5H + 96 f'H$6H = . (96) 



Solving (95) for 5H and plugging into (96) gives an equation for 
5$ + 3H 5<j) + 



V" - ZMK 2 f 2 H% + 24/"^ = Q 



h + 96K 2 f l2 H* 
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(97) 



where we used (88) to eliminate V . Specializing to our model we set 

V{<f>)= Pvac -J<t>, /(</>) = ^ 2 2 - 



(98) 



Then (97) reduces to (for p vac 3> p^., an d using (36) to eliminate the vev) 



4 

™ Pvac 



h u vi+4(e//i)«v 



5cj) = 0. 



(99) 



This means the low-curvature de Sitter solution is unstable. Assuming p vac ~ 
mp and ^/h < 1 the factor in parenthesis is 0(1). So the timescale for the 
instability is 



1/2 



Ho 1 . 



(100) 



Provided £/7t < 1 the solution will be metastable for more than a Hubble 
time, and Hubble friction will play an important role in the evolution of (f>. 
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